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Abstract 

We show that a nontrivial graph isomorphism problem of two undi- 
rected graphs, and more generally, the permutation similarity of two 
given n x n matrices, is equivalent to equalities of volumes of the in- 
duced three convex bounded polytopes intersected with a given sequence 
of balls, centered at the origin with radii t$ G (0, \/n — 1), where {ij} 
is an increasing sequence converging to yn — 1. These polytopes are 
characterized by n 2 inequalities in at most n 2 variables. The existence 
of fpras for computing volumes of convex bodies gives rise to a semi- 
fpras of order 0*(n 14 ) at most to find if given two undirected graphs 
are isomorphic. 
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1 Introduction 

Let Gi = (y,Ei),G 2 = (V,E 2 ) be two simple undirected graphs, where V is 
the set of vertices of cardinality n and E 1 ,E 2 C V x V the set of edges. Gi and 
G 2 are called isomorphic if there exists a bijection a : V — > V which induces 
the corresponding bijection a : E\ — > E 2 . The graph isomorphism problem, 
abbreviated here as GIP, is the computational complexity of determination 
if G\ and G 2 are isomorphic. Clearly the GIP in the class NP. It is one 
of a very small number of problems whose complexity is unknown [HI [TO] . 
For certain graphs it was known that the complexity of GIP is polynomial 

P El [3 S El [IE]. 
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The graph isomorphism problem is a special case of permutational sim- 
ilarity of two nxn real values matrices A, B G IR ?1X ' 1 , or more generally 
nxn matrices with entries in any ring with identity. Namely, let V n C M nxn 
be the group of permutation matrices. Does there exists P G V n such that 
B = PAP T 1 

Using the notion of coherent algebras, as a tool to identify nontrivial pairs of 
matrices A, B G M. nxn which may be permutationally similar [TJ, we first show 
that a nontrivial permutational similarity of A, B G M. nxn can be polynomially 
reduced to an isomorphism problem of two regular undirected connected multi- 
graphs with the same degree, (self- loops allowed), and the same characteristic 
polynomial. Assume that these two graphs represented by symmetric A, B G 
S n (Z + ) its rows sum equal to N. 

Denote by 1 = (1, ... , 1) T , and let M n C R nxn be the space of matrices 
with zero row and column sum: 

XI = X T 1 = 0, X = [ Xij ] G R nxn . (1.1) 

Note that M n is (n — l) 2 dimensional subspace of M. nxn . Denote by fl n cR" x " 
the convex set of doubly stochastic matrices. Note that Y G Q n if and only if 
Y = X + J n , J n := ^H T , where X G M n and each entry of X satisfies t > —-. 
It is easy to see that := V tr XX T < \Jn—\ if — J n < X G M n , and 

equality holds if and only if X = P — J n for some P G V n - 

For any S,T G W nxn we define the following subspace of matrices and a 
corresponding bounded polytope: 

P (5,T) := {X G M n , SX-XT = 0}, P(S,T) := {X G P (S,T), X > -J n }. 

(1.2) 

Denote by ||P(5', T)\\ F := maxxep(5,T) II^IIf, the radius of P(S,T). Thus two 
regular multi undirected graphs G, H, with the the same number of vertices 
and edges, are isomorphic if and only if ||P(A, B)\\ F = \/n — 1, where A, B are 
the representation matrices of G, H respectively. However, it is known that 
finding the radius of a convex set is iVP-hard [15]. (Note that the diameter of 
a balanced convex set K, i.e. — K = K, is twice its radius.) 
The main result of this paper is 

Theorem 1.1 Let G,H be two regular multi undirected graphs, with the 
the same number of vertices n and edges e. Denote by A, B G Z™ xra the 
representation matrices of G, H respectively. The the following statements 
are equivalent. 

1. G and H are isomorphic. 

2. The dimension of the convex sets P(A, A), P(A, B), P(S, B) are equal. 
Furthermore, for each t G (0, \Jn — 1] the volumes of the intersection of 
the above three polytopes with the ball of radius t centered at are the 
equal. 
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3. The dimension of the convex sets P(A, A),P(A, B),P(B, B) are equal. 
Furthermore, for a given sequence of balls ti G (0, \/n — 1), where {U} is 
an increasing sequence converging to y/n — 1, the volumes of the inter- 
section of the above three polytopes with each ball of radius ti centered at 
are the equal. 

The main argument of the proof of this is theorem follows straightforward 
from the observation that I n — J n G P(A, A), i.e. ||P(A A)\\ F = y/n — 1. 

We show that for the convex sets P(A, A), P(A, B), P(B, B) one can ap- 
ply the known results, which give fully randomized polynomial approximation 
scheme for computing the volumes of the intersection of these set with a ball 
of radius t G (K y/n — 1), e.g [U [151 El ESj. Combining these results we ob- 
tain some algorithms for testing the volume conditions given by Theorem 11.11 
Recall that the problem of finding the exact volume of a polytope in IR m , given 
by a polynomial number of affine inequalities in m, is #P-hard [5]. We hope 
that this approach will lead to a fpras to determine if given two graphs are 
isomorphic. 

We now summarize briefly the contents of this paper. In §2 we discuss 
the notion of coherent algebras and their relations to the graph isomorphism 
problem. In §3 we construct the polytopes P(A, A), P(A, B), P(B, B) which 
are intimately related to a permutational similarity of A, B, assumed to be 
scaled doubly stochastic with the same row sums. We give an outline of the 
proof of Theorem 11.11 In §4 we outline a semi-fpras to find if given A, B are 
permutationally similar, which is based on the fpras for computing volume of 
convex sets. 

2 Coherent algebras 

A subalgebra C C W nxn is called a coherent algebra, if it is closed under 
the transposition and entry-wise multiplication of two matrices, and contains 
I n = [<%] , J n = [^] , the identity matrix and the doubly stochastic matrix 
with equal entries. (Denote by A o B = [ay] o [fey] the entry-wise product 
[dijbij].) We now briefly survey the main properties of coherent algebra used 
her. Our main source is our paper [7|. Additional references for the properties 
of coherent algebras cited explicitly, where needed. 

A trivial coherent algebra is an algebra of dimension 2 spanned by /, J. Co- 
herent algebras of dimension 3 are induced either by strongly regular graphs, 
or by Hadamard matrices. A coherent algebra C has a canonical basis consist- 
ing of (0, 1) matrices E±, . . . , Ed- Each Ei is either symmetric or asymmetric, 
i.e. Ei o Ej = 0, and balanced, i.e. the nonzero rows and columns of Ei 
are equal to r« and q respectively. Furthermore, Ei o Ej = for i ^ j, and 
Y? i=l Ei = 11 T . C is characterized by the tensor T(C) = [t hj , k ] G Z d + Xdxd . 
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Any A G M nxn induces the minimal coherent subalgebra C(A) C M" x ™ 
containing A. One finds in polynomial time the canonical basis E\ , . . . , G 
{0, l} nxn of C(A). ([71 Lemma 3.1] yields that one needs at most Yin 10 flops.) 
If B G R nxri is permutationally similar to A, it follows that C(B) is strongly 
isomorphic to C(A). So C(B) has the canonical basis Fi, . . . , Fd, such that 
Fi = PEiP T , i — 1, . . . , d for a corresponding P eV. These equalities induces 
the strong isomorphism t : C(A) — > C(-B) given by t(-Ei) = i = 1, . . . , d. 

Thus for A, B G l"*" 1 to be permutationally similar we must have an 
isomorphism t : C(A) — > C(B), such that = Fj, z = 1, . . . , n. (4 is an 

isomorphism of two algebras, which preserve the transposition and entry-wise 
multiplication, i.e. i{U T ) = i(U) T ,i(U o V) = l(U) o l(V).) In particular, 
T(C(A)) = T(C(B)), and this condition is essentially equivalent to isomor- 
phism ofC(A) and C{B). 

The existence or nonexistence of such isomorphism is determined in a poly- 
nomial time in 0(n w ). It is possible that the isomorphism of coherent algebras 
does not imply the strong isomorphism. 

Theorem 2.1 Let A, B G M nxr \ Assume that the coherent algebras C(A) 
and C(B) are isomorphic, i.e. 1 : C(A) — ► (B) is an isomorphism of coherent 
algebras. Then there exists A\ G C(A), B\ G C(-B) wzt/i the following properties: 
i(Ai) = l(Bi); A\ and Bi are symmetric matrices with positive integer entries 
whose values are less than n 3 ; each row sum of A\ and Bi is equal to N ; A\ 
and B\ have the same characteristic polynomial. Furthermore, A and B are 
permutationally similar if and only if A\ and B\ are permutationally similar. 

Outline of Proof. Assume Ei, . . . , E^, h > 1 are all the diagonal ma- 
trices in the canonical basis E±, . . . , Ed- Consider A2 := Ylt=h+i m i^ii where 
m 2 , . . . , rrid G N satisfy the condition rrii 7^ rrij unless Ej = Ej. (If Ej = Ej we 
let rrti = rrij.) The number of distinct integers in {m^+i, . . . , rrid} is p < n ( n ~ 1 > . 
Hence we can assume that the set of distinct integers in {m h+ i, . . . ,m d } is 
{1, . . . ,p}. Let N — 1 be the maximal row sum of A 2 . Then A\ = A + D, 
where D is the diagonal matrix such that each row of A2 equals to N. The 
results of [7J imply that A\ G C(A). Set B\ = i{A\). Then all other claims of 
the theorem follow straightforward from the results in [7J. □ 

Note that A\ , B\ are representation matrices of two regular undirected 
multi-graphs G, H with self loops, with the same numbers of vertices, edges, 
and the same characteristic polynomials. Furthermore, G and H are con- 
nected. In the rest of the paper we assume that A — A%, B — B\. 

It is possible to show, that by increasing the entries of A%, which are still are 
of order 0(n K ), that in addition to the above conditions on Ai, A\ is generic in 
C(A). That is, the multiplicity of each eigenvalue of A% is the minimal possible 
for any symmetric matrix S G C(A). 

More generally, any set of Ax, . . . , Af. G IR nxn induces a minimal coherent 
subalgebra C(Ai, . . . , A^) C R nx " which contains these matrices. It is obtained 
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by the following process. First, express each matrix A4 as a unique linear 
combination of (0, 1) matrices with pairwise distinct coordinates. This gives 
rise to T u . . . ,T k G {0, l} nxn . (At the first step T\ = I n} T 2 = 11 T .) By 
considering the subspace spanned by all nonzero (0, 1) matrices of the form 
TioTj, T^oTj, Tj 7 °Tj , i, j = 1, . . . , k we obtain a subspace Ui C M nxn spanned 
by (0, 1) matrices R\, . . . ,Ri with disjoint support such that their sum is equal 
to 11 T . We now consider the set of matrices RiRj, i, j = 1, . . . , I, whose span 
U 2 includes Ri, . . . , R[. Apply the previous algorithm to these I 2 matrices to 
obtain a (0, 1) basis in U2 which is a refinement of the basis Ri, . . . , Ri. After 
p < n 2 steps we will have that U p = U p+ i. Then C(A 1 , . . . , A k ) = XJ P . 

We say that the set Ai, . . . ,A k G IR nxn is permutationally similar to 
B u ...,B k G R nxn if B { = PAiP T for i = 1, . . . , k and some P G V n . As 
in the case k = 1 the nontrivial permutational similarity induces an isomor- 
phism of the coherent algebras 1 : C(A\, . . . , A k ) — > C(Bi, . . . , B k ), such that 
i(Ai) = Bi,i = 1, . . . , k. As in the case k — 1, the problem of nontriv- 
ial permutation similarity of Ai, . . . , A k and B±, . . . , B k can be reduced to 
permutational similarity of two symmetric scaled doubly stochastic matrices 
A, B G M nxn with the same characteristic polynomial. 

3 Convex polytopes associate with GIP 

Let Q n C M" xri be the convex set ofnxn doubly stochastic matrices. Recall 
that tt n = {X e M„, X > -J n }. We say that A G M+ Xn is a scaled doubly 
stochastic matrix if A = aC for some a > and C doubly stochastic. 

Lemma 3.1 Let A G M™ xn ,n > 2 be an irreducible symmetric scaled dou- 
bly stochastic matrix. Assume that A has fi + 1 > 2 distinct eigenvalues 
A > Ai > . . . > A M , where rrii is the multiplicity of A, ; for i = 0, 
(mo = I.) Then the dimension of the subspace Po(A, A) and the polytope 
P(A,A), given by FTE . is A := YTi=\ m \- 

Assume that B G IR™ xn is an irreducible symmetric scaled doubly stochastic 
matrix having the same row sums as A. Then the following are equivalent 

1. A and B have the same characteristic polynomial. 

2. The three subspaces Po(A A), Po(A B), Pq(B, B) have the same dimen- 
sion. 

Proof. Since A is irreducible, its Perron- Frobenius root, Ao is the largest 
eigenvalue of multiplicity 1. As A is scaled doubly stochastic, A = X C for 
some symmetric stochastic matrix C. So AI = \ 1. Choose an orthonormal 
basis of K n consisting of orthonormal eigenvectors x = ^1, x 1; . . . , x n _ 1; 
corresponding to the eigenvalues Ao > Ai > . . . > A M . Let Q G M nxn be the 
orthogonal matrix whose columns are x , . . . , x n _i. Then Q T AQ is the block 
diagonal matrix ©f =0 Aj/ mr Observe that 

Q T M n Q = {0 lxl © Z, Z G (3.1) 
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Recall next that any commuting matrix with Q T AQ has the block diagonal 
form ®f =0 Ui where U G M m * xm * for i = 0, . . . , fi. Hence 

Q T P (A, A)Q = {0 lxl ©f =1 Ui G M m * xm ' , « = 1, . . . , fj}. (3.2) 

Thus dimP (A, A) = A. Since nxn G P (A, A) is an interior point of P(A, A) 
it follows that dimP(A, A) = A. 

Assume first that B has the same characteristic polynomial as A. Since 
A, B are symmetric, there exists an orthogonal matrix Qi, with he columns 
x ,yi, . . . ,y n -i, such that QjBQi = ©f =0 A;/ mi where t. Hence 

P (A, S) = P (A, A)Qj ; Q 2 = QiQ T , B = Q 2 AQ T 2 . (3.3) 

So dim P (A, 5) = dim P Q (B, 5) = A. 

Assume that dimP (A, A) = dimP (A, B) = dimP (-B, B). Hence the di- 
mension of the following three subspaces in M nxn : {X, AX— XA = 0}, {X, AX— 
XB = 0}, {X, BX — XB = 0} are equal. Therefore A and B are similar 

□ 

Proof of Theorem 11.11 Observe first that if U, V G M nxn are orthogonal 
matrices, then the transformation X i— > UXV is an orthogonal transforma- 
tion on M nxn . In particular, the ball B(0, t) C M raxn , centered at of radius t 
satisfies the equality UB(0, i)V = B(0, t) for any orthogonal U, V. 

Assume first that A, B are similar, i.e B = Q 2 AQ 2 Y for some orthogonal 
Q 2 . Then for any t > 

B(o, t) n P (A, B) = (B(o, t) n P (A, A))Ql, 

B(0, t) n P (B, 5) = Q 2 (B(0, t) n P (A, A))Qj. 

Suppose furthermore that A and I? are permutationally permutationally sim- 
ilar, i.e Q 2 = P G P n . (So M holds.) Then ([313]) yields that P(A,B) = 
P(A, A)P T , P(B, B) = PP(A, A)P T . In particular 

B(0, *) n P(A, B) = (B(0, t) n P(A, A))P T , 
B(0, t) n P(B, 5) = P(B(0, t) n P(A, A))P T . 

Hence, above intersections have the same volume for any t > 0. This proves 
the conditions [2J|3J 

Recall that ||P(A, B)\\ F < \Jn — 1, and equality holds if and only if ("P n — 
J n )flP(A, P) 7^ 0, i.e. A and B are permutationally similar. Observe next that 
si nce I n - J n G P(A, A)(~)P(B, B) it follows that ||P(A,A)|| F = \P(B,B)\\ F = 
yjn - 1. Hence, the volumes of B(0, t) (1 P(A, A), B(0, t) H P(B, B) increase in 
the interval (0, \Jn — 1). 

Assume that [3] holds. Since the volumes of B(0, tA H P(A, A), i = 1, . . . , 
form an increasing, it follows that the volumes of B(0, tA fl P(A, B), i — 1, . . . , 
form an increasing sequence. Hence ||P(A, B)\\ = y/n — 1. So A and B are 
permutationally similar. Clearly, the same arguments imply that the condition 
|2] implies the permutational similarity of A and B. □ 
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4 A semi-fpras for graph isomorphism 



We now point out how to apply the existing fully polynomial randomized ap- 
proximation schemes to compute a volume of a convex sets, e.g. [TT] . To 
do that it would be convenient to map the three convex polytopes P(A,A), 
P(A, B), P(B, B) of dimension A, to one ambient space M A by by three dif- 
ferent linear transformations: 

Tt : P (A A) - R A , T 2 : P (A, B) - R A , T 3 : P (A, A) - M A (4.1) 

such that each preserves the inner product. We demonstrate for T 2 . Choose 
an orthonormal basis W\, . . . , Wa in Po(A -B). Then T(Wi) = (5u, . . . , 5 n A) T , i = 
1, . . . , A. It is straightforward to show that 7i(P(A, A)),T 2 (P(A, B)),T 3 (P(B, B)) 
are polytopes Xi,X2,X3 C M AxA , which are given as follows. There exists 
nonzero vectors u^iy, . . . , VL( n ,n),i G K A , I — 1, 2, 3 such that 

X = {x G M A , uj^x >~i,j = l n}, Z = 1, 2, 3. (4.2) 

One can compute the vectors uu,j),i in polynomial time in n. Note that the 
inequality > — ^ is equivalent to uX^ z x > — - in the orthonormal basis of 
the corresponding linear space Po(^4, A), Po(^4, B), Pq(B, B). 

As we explain in the next section the permutational similarity of A and B 
is equivalent to the existence of an orthonormal matrix O G M AxA such that 
0{u(i,i),i, . . . , U( n)n)i i} = {u (li i )i2 , . . . , U(„ )n ) j2 }. (Similar statement holds for 

{ u (l,l),3 ; • • • i u (n,n),3}, {U(l,l),2, ■ • ■ j U( n)n ) )2 }.) 

It is trivial to see that the polytopes P(A,A),P(A,B),P(B,B) contain 
the ball of radius - centered at the origin. Theorem 11.11 yields that A and B 
are permutationally similar if 



vol(B(0. i)Xi) = vol(B(0, t)X 2 ) = vol(B(0, t)X 3 ) for each t G (-, y/n- 1). 

n 

(4.3) 

We now suggest a probabilistic test of the above equalities, for a finite 
number of values t G (-, y/n — 1 with a relative error e and with probability 
1 — rj. In the random algorithms suggested in [U [151 adopted to find the 
volumes of Xi, i — 1,2,3, one considers the intersection of the sequence of balls 
of radii: 

2^ , 

tj = —, j = 0,...,M= \Nlog 2 nVn^T\. (4.4) 

Here can be chosen as A, as in [U [151 [11], or if we want more points we can 
take N = n c A for some c> 0. Let X i(< := B(^) n X i5 j = 0, . . . , M, i = 1, 2, 3. 
For each X^ one generates p = 400e: _2 A^log random points from certain 
distribution, e.g. [HI §6]. Then v ^ x,, j|j is estimated by the fraction of 
number of the sampled points in X^ to the number of the sampled in X^j 
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which are in Xj_! j. By Theorem II .![ if A and B are permutationally similar, 
we must have the equalities 



vo\(X jtl ) vol(X ii2 ) vol(X i)3 ) 



vo^X^i) = vol(Xj j2 ) = vol(Xj j3 



voKX^x,!) vol(X,_ li2 ) vo\(X 3 ^ 3 ) 

(4-5) 

for j — 1, . . . , M. Hence, in our process of estimating the volumes of Xi, X 2 , X3 
we test the above equalities within relative error e. If all the above equalities 
hold within the relative error e, we declare that the matrices A, B are e, rj 
permutationally similar. If one of the equalities in (14. 5p fails with respect to 
relative error e, we have two options. Either declare that the matrices A, B 
are not e, r) permutationally similar, or retest this equality with a smaller e 
and rj. If all the retested equalities hold, then we declare that A, B are e, 77 
permutationally similar. Otherwise, we declare that the matrices A, B are not 
e, 77 permutationally similar. 

Note that each oracle query if a point x 6 B(i) fl X; needs An 2 < n 4 flops, 
since the dot product in M A need A flops. Since the randomized algorithm 
suggested in [TTJ is of order 0*(A 5 ) we see that the our algorithm for checking 
the graph isomorphism, or permutational similarity of A, B G lR nx?1 , is of order 
0*(n 14 j at most. 
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